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ABSTRACT 

The  aim  of  this  dissertation  is  twofold.  The  early  chapters 
consist  of  an  expository  review  of  the  historical  development  of 
relativistic  cosmology  since  Einstein  suggested  the  general  theory  of 
relativity  in  1916.  One  chapter  is  devoted  to  developing  the 
corresponding  Newtonian  theory  and  it  is  seen  that  many  results  are 
analogous  with  the  relativistic  case.  In  the  latter  part  of  the 
thesis  some  new  solutions  of  the  cosmological  field  equations  are 
suggested;  one  of  these  solutions  is  shown  to  characterize  an  empty 
rotating  universe  in  which  test  particles  (nebulae)  of  negligible 
mass  recede  from  the  axis  of  rotation.  Hence  this  model  combines 
features  of  the  de  Sitter  and  Godel  type  universes. 


■ 


ACKNOWLEDGMENTS 


I  would  like  to  express  my  appreciation  to  Professor 
J.  R.  Trollope  for  suggesting  the  problem  in  the  latter  part  of 
this  dissertation  and  for  many  helpful  and  illuminating  discussions. 
I  also  wish  to  thank  my  wife,  Laurel,  for  her  diligence  in  preparing 


the  manuscript. 


# 


TABLE  OF  CONTENTS 


Page 

ABSTRACT  .  i 

ACKNOWLEDGEMENTS  .  .........  ii 

CHAPTER  I  INTRODUCTION  .  1 

CHAPTER  II  NEWTONIAN  COSMOLOGY  .  6 

CHAPTER  III  RELATIVISTIC  COSMOLOGY:  ROBERT SON- WALKER 

METRIC . 15 

CHAPTER  IV  NON-STATIC  UNIVERSES  .  26 

CHAPTER  V  THE  GODEL  SOLUTION:  A  ROTATING  UNIVERSE.  36 

CHAPTER  VI  SOME  NEW  SOLUTIONS  OF  THE  COSMOLOGICAL 

FIELD  EQUATIONS .  42 

CHAPTER  VII  EXAMPLES  AND  DISCUSSIONS  OF  RESULTS  ...  55 

BIBLIOGRAPHY  . . .  63 


* 

- 


CHAPTER  1 


Introduction 

Astronomical  investigations  of  the  entire  observable  universe 
have  led  to  two  results  which  are  of  primary  importance  to  modern  cosmo¬ 
logical  theories o  They  are  the  uniformity  of  the  large  scale  distribution 
of  extragalactic  nebulae  about  our  galaxy  and  the  observation  of  a  red- 
shift  in  the  spectral  lines  of  light  from  the  nebulae. 

While  individual  nebulae  are  to  be  found  singly,  in  pairs  or 

in  clusters  yet  on  a  very  large  scale,  say  that  of  clusters  of  galaxies, 

the  distribution  is  quite  uniform.  The  observable  region  is  both 

isotropic  and  homogeneous  i.e.  it  is  much  the  same  everywhere  and  in 

all  directions.  Hubble  [1]  notes  that  on  the  average  the  nebulae  are 

approximately  two  million  light  years  apart  i.e.  the  interval  between 

the  nebulae  is  about  two  hundred  times  the  mean  diameters.  The  mean 

-31 

density  of  ponderable  matter  is  calculated  to  lie  between  10  and 
-27  3 

10  g/cm  .  Of  course  this  estimate  of  mass  density  neglects  the 
presence  of  intergalactic  dust  which  could  conceivably  account  for  much 
of  the  matter  in  the  universe. 

The  displacement  toward  the  red  of  spectral  lines  in  the  light 
received  from  the  nebulae  is  known  to  increase  with  the  apparent  faintness 
of  the  source.  Since  astronomical  distances  are  determined  by  apparent 
faintness  it  is  evident  that  red-shifts  increase  with  the  distances  of  the 
observed  nebulae,  the  relation  being  approximately  linear.  The  most 
popular  physical  interpretation  of  red-shift  is  to  regard  it  as  a  Doppler 
effect,  in  which  case  it  follows  that  the  nebulae  are  receding  from  our 


. 

. 


2 


stellar  system  with  velocities  that  increase  directly  with  distance. 

This  is  the  origin  of  the  "expanding  universe"  theory.  It  is  worth 
noting  that  there  are  other  mathematically  consistent  interpretations 
of  nebular  red-shifts,  for  example  Sen  [2]  shows  it  to  be  a  consequence 
of  an  inherent  geometrical  property  of  a  static  universe  whereas 
Gillieson  [3]  explains  it  as  a  result  of  the  inhomogeneities  in  the 
gravitational  field. 

The  Cosmological  Principle  is  an  extension  of  the  experimental 
conclusion  that  the  observable  universe  is  uniform.  The  principle  may 
vary  considerably  from  one  theory  to  another  but  essentially  it  remains 
a  fundamental  postulate  which  asserts  the  uniformity  of  the  entire 
universe;  such  an  hypothesis  facilitates  the  construction  of  the  many 
logically  reasonable  world  models.  As  a  consequence  of  the  cosmological 
principle  we  conclude  that  our  stellar  system  cannot  be  regarded  as 
occupying  a  privileged  position  in  the  universe;  the  assumption  of 
homogeneity  implies  that  any  two  positions  in  the  universe  are  equivalent 
for  the  purpose  of  astronomical  observation  and  physical  experimentation. 
Isotropy  further  implies  that  there  is  no  favoured  direction  in  space  and 
that  regardless  of  the  direction  in  which  the  telescope  is  trained  the 
observer  will  find  roughly  the  same  matter  distribution.  Hence  the 
cosmological  principle  can  be  regarded  as  an  extrapolation  to  the  universe 
as  a  whole  of  the  observations  of  an  observer  in  a  limited  region  of 
space-time. 


As  a  starting  point  for  specific  theories  it  is  necessary  to 


formulate  the  principle  in  more  precise  mathematical  terms.  The 


' 
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Robertson-Walker  line  element  [4]  describes  the  geometry  of  space-time 
in  a  large  class  of  relativistic  universes.  One  form  of  the  cosmological 
principle  which  leads  to  its  derivation  involves  three  fundamental 
hypotheses : 

1.  Viewed  on  a  cosmological  scale  space-time  can  be  separated 
into  a  "cosmic  time"  t  and  a  set  of  privileged  space-like 
hypersurfaces,  often  referred  to  as  the  substratum. 
Corresponding  to  each  value  of  t  there  is  one  and  only  one 
such  hypersurface.  The  "cosmic"  or  "global"  time  t  is 
actually  the  proper  time  measured  from  a  given  hypersurface 
along  the  timelike  geodesics  orthogonal  to  it. 

2.  The  hypersurfaces  t  =  constant  are  isotropic. 

3.  Any  two  points  on  a  hypersurface  t  =  constant  are  equivalent 

The  Steady  State  theory  of  Bondi,  Hoyle  and  others  [5]  adds  the 
further  restriction  that  the  geometry  of  the  three-dimensional  space-like 
hypersurfaces  must  be  time  independent.  This  assumption,  however,  leads 
to  the  contradiction  of  the  principle  of  conservation  of  matter,  a  result 
which  renders  the  theory  unacceptable  to  a  considerable  portion  of  the 
scientific  world. 

Early  attempts  in  the  nineteenth  century  to  solve  the 
cosmological  problem  within  the  framework  of  Newtonian  theory  met  with 
little  success,  the  failure  being  due,  for  the  most  part,  to  the  ad  hoc 
assumption  that  the  matter  distribution  in  the  infinite  Newtonian  universe 
is  not  only  uniform  but  also  at  rest  in  a  given  inertial  frame  (the  ether) 


. 
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However  this  assumption  gives  rise  to  Olber^s  paradox  {6]  according  to 
which,  in  this  universe,  there  is  no  reason  for  the  sky  to  be  dark  at 
night!  This,  and  other  difficulties  inherent  in  the  above  model 
discouraged  further  interest  in  cosmology  until  Einstein  started 
applying  his  General  Theory  of  Relativity  to  the  problem.  In  1934 
Milne  and  McCrea  [7]  reopened  the  question  of  Newtonian  cosmology  and 
by  suitable  modification  of  the  Einstein  field  equations  were  able  to 
develop  a  theory  which  is  in  many  respects  analogous  to  the  relativistic 
theory.  One  of  the  obvious  benefits  of  such  a  Newtonian  conception  of 
the  universe  is  that  the  various  terms  are  more  readily  given  a  physical 
interpretation  which  may  then  pass  over  into  the  relativistic  case  by 
analogy.  It  may  be  pointed  out,  however,  that  while  this  is  a  very 
interesting  mathematical  exercise,  some  of  the  tenets  of  Newtonian 
dynamics  have  been  clearly  proven  to  be  untenable  in  the  light  of  later 
work,  particularly  the  special  theory  of  relativity;  hence  in  order  to 
obtain  any  reasonable  hypothesis  about  the  state  of  the  universe  we  must 
turn  to  some  post-Newtonian  theory,  general  relativity  being  the  most 
widely  acceptable  to  date. 


The  Einstein  field  equations  are 


R..  -  ^  g.  .  R  +  Ag..  =-  8ttkT.  . 

ij  2  6ij  6ij  ij 


(1.1) 


where  g..  is  the  metric  tensor,  R. .  is  the  Ricci  tensor,  R  =  g1^  R. . 

ij  iJ 


A  and  k  are  universal  constants  and  T. .  is  the  energy-momentum  tensor. 

IJ 

The  Newtonian  limit  analog  of  (1.1)  is 


V^<J>  +  Ac^  =  4ttkp  , 


(1.2) 


_  i  s rj  -  .  .  .  ■ ;  mnyb 
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where  p  is  the  mass  density  and  <f>  is  gravitational  potential. 


Einstein  originally  introduced  the  cosmological  constant  A 

in  order  to  obtain  a  static  cosmological  solution.  In  1922  Friedmann 

found  some  non-static  solutions  of  the  Einstein  equations  with  A  =  0 

which  satisfied  the  cosmological  principle.  In  view  of  this  and  other 

considerations  Einstein  eventually  abandoned  the  cosmological  term  in 

the  field  equations.  There  is  considerable  difference  of  opinion  about 

the  value  of  including  A  in  modern  cosmological  theory.  However  special 

•  • 

relativity  requires  that  T  J  ^  =  0  and  since  it  can  be  shown  that  the 
quantity 


Rlj  -  |  glj  R  +  Aglj 


is  the  most  general  second-rank  tensor  with  zero  divergence  which  can  be 
constructed  entirely  from  the  metric  tensor  and  its  first  and  second 
derivatives,  it  seems  that  the  inclusion  of  A  is  desirable.  Also  once 
solutions  of  (1.1)  are  found  then  by  considering  A  as  a  variable  parameter 
we  can  admit  a  wider  range  of  model  universes.  It  follows  that  A  may  be 
valuable  in  explaining  hitherto  unsuspected  small  effects  when  more 
sophisticated  observational  techniques  are  developed  and  more  detailed 
data  become  available.  It  is  worth  noting  that  in  "local"  celestial 
mechanics,  for  example  on  the  scale  of  our  planetary  system,  the  effects 
of  A  are  so  small  as  to  be  experimentally  undetectable.  The  field  equations 
may  then,  without  loss  of  generality,  be  written  as 


R.  . 
iJ 


1 

2 


R  =  -  8ttkT  .  . 

ij 


(1.3) 


However  on  the  cosmological  scale  the  effects  of  even  a  very  small  A 
could  be  of  major  significance. 


T;'\v  ■c.i. 
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CHAPTER  2 


Newtonian  Cosmology 

In  view  of  our  previous  remarks  it  is  necessary  to  investigate 
solutions  of  the  field  equation 

V^cp  +  Ac^  =  4ttkp  o  (2.1) 

In  this  theory  there  is  of  course  an  even-flowing  Newtonian  time  t  which 
is  the  same  for  all  observers.  The  appropriate  form  of  the  cosmological 
principle  is  that  the  three-dimensional  hypersurfaces,  called  the 
substratum,  corresponding  to  each  fixed  value  of  t,  are  composed  of  a 
homogeneous  fluid.  Any  two  observers  moving  with  the  substratum  will 
obtain  the  same  results  provided  they  conduct  their  experiments  at  the 
same  time. 


Consider  now  two  typical  observers  0  and  O’  ,  moving  with  the 
substratum,  both  of  whom  are  studying  the  motion  of  a  general  point  P 
of  the  fluid.  Let  r  and  r'  be  the  position  vectors  of  P  with  respect  to 
0  and  0  respectively.  Let  the  vector  00  be  denoted  by  a. 


o 
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Let  v(r,t)  be  the  velocity  of  P  as  observed  by  0  and  v'(r’,t)  be  the 
velocity  of  P  as  observed  by  O' .  It  follows  from  simple  vector  addition 
that 

v(r,t)  =  v(a,t)  +  v* (r ' ,t)  .  (2.2) 

Now  suppose  that  P*  is  another  point  of  the  substratum  whose  position  with 
respect  to  O’  is  similar  to  the  position  of  P  with  respect  to  0;  then  by 
the  cosmological  principle  the  velocity  of  P  as  measured  by  0  must  equal 
the  velocity  of  P'  as  measured  by  O'.  Thus  we  get 

v(r, t)  =  v* (r,t)  ,  (2.3) 

which  must  hold  for  any  r.  Combining  this  equation  with  (2.2)  gives  rise 
to  the  result 

v(r-a,t)  =  v(r,t)  -  v(a,t)  .  (2.4) 

This  last  equation  shows  that  v  is  a  linear  function  of  r,  hence,  in  terms 
of  components  we  may  write 

v1  =  A1.^  ,  (2.5) 

3 

where  A1^  is  a  tensor.  However,  the  characteristic  equation  of 

A1 .r .  =  Ar.  ,  (2.6) 

J  i  J 

is  a  cubic  and  therefore  has  at  least  one  real  root.  Hence,  there  is  a 
direction  r  such  that  v  is  directly  away  from  the  observer.  But  since  the 
space  is  isotropic,  this  must  be  true  for  all  directions  and  we  conclude 
that 

v  =  A(t)  r 


(2.7) 


. 


'  .  1  '  >-*>*<?  8 
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To  solve  (2.7)  define  a  new  function  R(t)  by  the  relation 


A(t)  .  1  dR  =  R 
M  '  R  dt  ‘  R 


(2.8) 


Integration  then  yields 


r(t)  =  R(t)r 


0 


(2.9) 


where  r, 


0 


r(tQ)  and  R(tQ)  =  1. 


In  view  of  our  assumption  that  the  substratum  is  composed  of 
a  perfect  fluid  we  are  now  in  a  position  to  introduce  two  important 
equations  viz.  the  equation  of  continuity 


ia 

3t 


+  V* (pv)  =  0 


(2.10) 


and  the  equations  of  motion  of  a  perfect  fluid  without  pressure 


dv 

dt  =-V+ 


(2.11) 


Note  that  in  accordance  with  the  cosmological  principle  p  is  a  function 
of  t  only  i.e.  p  =  p(t). 

Combining  equations  (2.7),  (2.8)  and  (2.10)  we  obtain 


lt+pfv,?=o  > 


(2.12) 


which  can  be  written  in  the  alternative  form 


dp  „dR 
P  ^ 


0 


(2.13) 


rin 
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Integration  yields  pR  =  constant  and  from  initial  conditions  we  conclude 


p(t)  = 


[R(t)] 


(2.14) 


Equations  (2.1)  and  (2.11)  lead  to  the  result 


oR  ,  »  2  .  1 

-  3g  +  Ac  =  4«pq-  -3 


(2.15) 


Multiplying  (2.15)  by  the  integrating  factor  RR  ,  integrating  and 
rearranging  we  get 


8tt<p 


R  - 


0 


1 

R 


Aci  2 

3  R 


"  B  , 


(2.16) 


where  8  is  a  constant  of  integration.  Now  put 


V (R)  =  - 


8tt<p 


0 


1 

R 


Aci  2 

3  R 


(2.17) 


to  get 

R2  +  V (R)  =  6  ,  (2.18) 

which  can  be  interpreted  as  the  energy  integral  if  R  is  considered  to  be 
the  coordinate  of  a  particle  moving  in  one  dimension  and  V(R)  is  identified 
with  potential  energy. 

Turning  our  attention  once  again  to  equation  (2.14)  it  can  be 

seen  that  with  our  present  interpretation  of  R  as  a  radial  coordinate  it 

follows  that  if  all  linear  dimensions  are  increased  by  a  factor  R(t)  then 

the  corresponding  density  decrease  must  be  proportional  to  ;  hence 

RJ 

this  equation  corroborates  a  result  which  is  intuitively  obvious. 


. 


<i  >  /  ;  .’rv  I  >1  q  m  70  i<  3  i  b 
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Putting  A  =  0  in  equation  (2.15)  we  obtain 

3R  ,  1 

-  —  =  4«P()  -3  . 

K. 


(2.19) 


Now  for  a  static  solution  we  require  R(t)  =  R(tg)  =  1  so  that  (2.19) 
becomes  simply  p  =  0  .  It  is  now  evident  why  the  assumption  that  there 
were  no  large  scale  motions  of  matter  in  the  universe  led  nineteenth 
century  cosmologists  to  conclude  that  it  was  necessary  to  consider  some 
generalization  of  Poisson's  equation  if  they  were  to  obtain  a  non¬ 
trivial  universe  with  a  positive  density  of  matter. 


We  are  now  in  a  position  to  examine  the  various  model  universes 
which  can  be  constructed  by  assigning  appropriate  values  to  the  arbitrary 
parameters  A  and  3.  The  simplest  case  is  A  =  3=0  when  the  equations  of 
motion  (2.16)  become 


8ttkp 
dR=  (— 


0)1/2 


dt 


(2.20) 


Integrating  and  choosing  R  =  0  when  t  =  0  we  get 


R  -  (6tticp0)1/3  t2/3 


(2.21) 


This  represents  a  universe  which  expands  monotonically  from  a  singular 
state  in  the  finite  past. 

As  in  relativistic  cosmology  the  admissible  models  range  from 
the  static  Einstein  universe  for  which  R  =  R^  to  the  oscillating  models 
in  which  R  =  0  at  t  =  0  after  which  there  is  a  period  of  expansion  up  to 
a  maximum  radius  R  =  Rq  followed  by  a  contraction  until  once  again  R  =  0  . 


(5  »  0  *lqoi±a  s?«c:o osd 

tJ*T  jo  xJicosb  ovr± 3lao<  i»  /j:Uv  -asovlfli.’  IftlvJtiJ' 

j-fl  X  O'  ,’OJCJ.  cCo-  A  n.t  *C  0  ftT*  .  1 
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This  important  class  of  cosmological  solutions  corresponds  to  the 
following  cases:  (i)  A  <  0,  all  3  ,  (ii)  A  =  0,  3  <  0, 

(iii)  A  >  0,  3  <  3q  where  3  =  3q  =  V(Rq)  when  R  =  R^  .  For  a 
detailed  treatment  of  the  different  models  see  Bondi  [6]  and 
Trautman.  [8], 


Finally  it  is  necessary  to  investigate  the  behaviour  of 


light  rays  in  a  Newtonian  universe.  One  must,  in  the  analysis, 
accept  the  pre-relativistic  concept  that  light  obeys  the  classical 
law  for  the  addition  of  velocities.  In  view  of  the  special  theory 
of  relativity  we  consider  this  to  be  a  false  premise;  nonetheless 
we  here  employ  it  for  the  sake  of  determining  the  properties  of  light 
in  the  universe  under  consideration.  The  results  indicate  that 
Newtonian  theory  gives  a  good  approximation  to  relativistic  cosmology 
as  long  as  local  velocities  are  small  compared  with  the  velocity  of 
light . 


In  accordance  with  the  cosmological  principle  the  local 


velocity  of  light  must  be  the  same  for  all  observers  at  a  given  time, 
though  it  may  be  a  function  of  time.  Applying  the  Newtonian  law  for 
the  addition  of  velocity  to  an  incoming  radial  light  ray  we  get 


(2.22) 


where  c(t)  is  the  local  velocity  of  light.  This  can  be  written  in  the 
form 


(2.23) 


e.  1;  sv  ffolJibbB  orfJ  it  v#al 
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We  can  regard  equation  (2.23)  as  expressing  the  law  of  propagation  of 
light  rays  in  a  Newtonian  universe.  In  order  for  this  law  to  be 
meaningful  within  the  framework  of  the  present  theory  it  is  necessary 
to  postulate  the  existence  of  an  ether,  expanding  with  the  universe, 
and  with  respect  to  which  "absolute"  velocities  are  measured.  For 
simplicity,  and  in  order  to  depart  from  the  relativistic  theory  as 
little  as  possible,  we  make  the  ad  hoc  assumption  that  c(t)  is  a 
constant.  Integration  of  (2.23)  then  yields 


where  r^  = 
of  a  light 
origin  r  = 


R, 


=  -c 


dt 

R(t) 


(2.24) 


r(t1),  r2  =  r(t2),  R1  =  R(t1) ,  R2  =  R(t2).  In  the  case 
ray  travelling  from  a  source  at  r(t)  to  an  observer  at  the 
0,  equation  (2.24)  becomes 


(2.25) 


Equation  (2.25)  determines  the  time  of  reception  t2  in  terms  of  the  time 
of  emission  t^  and  the  distance  of  the  source  at  the  time  of  emission 
i.e.  r^  .  Suppose  that  the  source  is  moving  with  the  substratum  and  that 
it  flashes  light  signals  at  successive  times  t^  and  t^  +  At-^  ;  suppose 
further  that  these  light  rays  are  received  at  the  origin  at  times  t2  and 
t 2  +  At2  respectively.  Since  the  source  is  moving  with  the  substratum 
the  left-hand-side  of  equation  (2.25)  viz.  r  /R^  is  independent  of  t^  and 


we  can  thus  write 


. 


(3)*’ 


13 


c 


t2+At2 


tl+Atl 


dt 

R(t) 


(2.26) 


Combining  equations  (2.25)  and  (2.26)  and  using  some  results  of  elementary 
calculus  we  obtain 


Atx  At2 

R(tx)  =  R(t2) 


(2.27) 


or 


At 


2 


At 


1 


(2.28) 


Now  suppose  At^  is  the  period  of  the  emitted  light  wave  and  At2  is  the 
period  of  the  received  wave.  Equation  (2.28)  immediately  suggests  that 
there  is  a  red-shift  only  if  >  R-^  ,  i.e.  if  the  universe  is  expanding. 
Writing  (2.28)  in  terms  of  wavelength  A(t)  we  get 


(2.29) 


Where  is  the  wavelength  of  received  wave  and  A^  is  wavelength  of 
emitted  wave.  The  shift  in  wavelength  z  is  defined  by 


z 


AA 

A 


(2.30) 


It  will  be  shown  later  that  this  is  identical  with  the  relativistic 
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expression  for  wavelength  shift.  It  is  interesting  to  note  that  Olber’s 
paradox,  referred  to  above,  can  be  resolved  by  considering  the  effects 
of  red-shift  in  an  expanding  universe. 
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CHAPTER  111 

Relativistic  Cosmology  :  Robert son-Walker  Metric. 

With  the  advent  of  the  general  theory  of  relativity  and  the 
early  successes  of  Einstein,  de  Sitter  and  others  in  constructing 
model  universes  within  the  framework  of  the  theory,  renewed  interest 
was  focused  on  the  cosmological  problem.  Indeed  modern  cosmology 
traces  its  origin  to  Einstein's  paper  of  1917  in  which  he  presented 
his  original  static  universe  [9].  As  previously  stated  most  of  the 
results  of  the  preceding  chapter  were  first  discovered  in  relativistic 
cosmology  theory,, 

The  structure  of  the  four-dimensional  space-time  continuum 
in  general  relativity  is  assumed  to  be  that  of  a  Riemannian  geometry. 
The  covariant  expression  for  interval  is 

ds^  =  g . ,  dx^  dx^  ,  (3.1) 

-^-3 

where  the  g„„'s  are  ten  functions  of  the  four  coordinates.  The  line 
ij 

element  ds  ,  which  gives  the  distance  between  two  neighbouring  points 
x1  and  x1  +  dx1  in  the  four  space,  is  assumed  to  be  an  invariant.  A 
further  condition  required  by  the  relativistic  theory  is  that  in  a 
sufficiently  small  region  of  space-time  there  exist  coordinates  such 
that 

2  2  2  2  2  2 
ds  =  -dx  -  dy  -  dz  +  c  dt  (3.2) 

i.e.  special  relativity  is  locally  valid.  The  structure  of  space-time 
is  determined  by  the  g „  which  in  turn  is  determined  by  the  field 


n  -  .  i  Ui  ig.t-ic  s 
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equations . 


The  field  equations  (1.3)  are  appropriate  for  a  cosmological 
model  filled  with  matter,  idealized  as  the  cosmological  fluid.  For 
convenience  we  may  choose  units  such  that  c  =  k  =  1  where  c  is  the 
velocity  of  light  and  k  is  the  universal  constant  in  (1.3).  We  then  get 


R.  . 
iJ 


-2gU  R 


8ttT  .  . 
iJ 


(3.3) 


Unless  otherwise  stated  we  shall  consistently  employ  these  units  in 
future.  Now  the  Einstein  tensor  R. .  -  ^  g..  R  is  constructed  entirely 
from  the  metric  tensor  and  its  derivatives  and  we  may  thus  regard  it  as 
specifying  the  geometry  of  space-time.  On  the  other  hand  the  energy- 
momentum  tensor  specifies  the  density,  momentum,  energy  and 

pressure  of  matter  in  the  universe.  Hence  equations  (3.3)  show  an 
explicit  relation  between  geometry  and  matter  content  and  in  this 
sense  the  Einstein  field  equations  are  an  embodiment  of  Mach's  principle. 

In  1883  Mach  [10]  postulated  that  the  inertia  of  a  body  is 
due  to  the  presence  of  all  other  bodies  in  the  universe.  In  1918 
Einstein  generalized  this  idea  and  put  forward  the  hypothesis  that  the 
geometry  of  space  is  determined  by  the  distribution  of  matter  in  the 
universe.  He  called  this  Mach's  Principle  and  was  of  the  opinion  that 
it  should  be  an  integral  part  of  general  relativity.  Hence  he  required 
that  the  T  should  uniquely  determine  the  g„  .  This  is  possible  only 
if  the  differential  equations  (3.3)  are  augmented  with  suitable  boundary 
conditions  at  infinity.  Einstein  demonstrated  that  no  such  boundary 
conditions  exist  and  in  an  attempt  to  overcome  this  difficulty  he 


. 

. 
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modified  equations  (3c 3)  by  introducing  the  cosmological  term  to  obtain 
the  new  equations 


R.  . 

ij 


o  g..  R  +  Ag. .  =  -8uT. . 
2  ij  ij  ij 


(3.4) 


For  A  >  0  equations  (3.4)  have  a  solution  for  which  mass 
density  p  is  constant  and  the  curvature  of  space  is  such  that  though 
unbounded  the  universe  is  finite.  This  is  the  Einstein  static  universe 
and  since  there  is  no  infinity  in  the  model  the  problem  of  finding 
boundary  conditions  does  not  arise.  Einstein  further  believed, 
erroneously,  that  there  existed  no  solutions  of  equations  (3.4)  with 
T  =  0  and  A  >  0  ;  this  would  imply  that  for  empty  space  there  can  be 
no  inertial  field,  an  argument  which  supports  Mach's  principle.  However, 
in  1917  de  Sitter  found  such  a  solution;  the  de  Sitter  universe  is 
empty  but  it  has  the  interesting  property  that  particles  of 
negligible  mass  can  be  shown  to  recede  from  each  other  with  steadily 
increasing  velocity,  thus  giving  rise  to  a  Doppler  effect.  The 
discovery  of  the  de  Sitter  model  was  one  of  the  factors  which  persuaded 
Einstein  to  abandon  the  cosmological  term.  One  is  forced  to  the 
conclusion  that  Mach's  principle  is  not  in  fact  an  intrinsic  property 
of  general  relativity  but  is  an  hypothesis  which  finds  one  form  of 
expression  in  the  relativistic  field  equations.  The  discovery  of 
Godel's  rotating  universe  [11]  in  1949  supports  this  conclusion  since 
it  demonstrates  that  by  using  the  same  energy-momentum  tensor  as  is 
assumed  in  the  Einstein  universe  one  can  find  a  totally  different 
geometrical  structure.  It  thus  appears  that  a  suitable  boundary  or 
global  condition  must  be  imposed  if  Mach's  principle  is  to  be 


» 
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considered  an  integral  part  of  general  relativity. 

In  relativistic  cosmology  the  nature  of  the  substratum  is 
specified  by  Weyl's  postulate  which  states  that  "the  particles  of  the 
substratum  (the  nebulae)  lie  on  a  bundle  of  geodesics  diverging  from  a 
common  point  in  the  past."  This  hypothesis  is  verified  by  observation 
to  be  substantially  correct  in  that  deviations  from  the  assumed  geodesic 
motion  are  found  to  be  random  and  of  small  velocity  compared  with  that 
of  light.  However  the  relative  velocities  of  the  nebulae  along  the 
geodesic  paths  may  be  very  large  and  will  thus  give  rise  to  the  red- 
shifts  observed  in  the  wavelength  of  light  received  from  the  nebulae. 

The  first  relativistic  world-models  assumed  the  cosmological 
principle  in  the  form  that  the  universe  is  both  homogeneous  and  has  a 
time- independent  geometrical  structure.  Clearly  space  may  be  regarded 
as  exhibiting  spherical  symmetry  about  any  observer  so  that  the  natural 
choice  of  line  element  is  the  sphericaily-symmetric  static  form 

,  2  A  ,  2  2J  2  2  .  2  2  ,  vJ  2  ,,  cN 

ds  =  -  e  dr  -  r  d0  -  r  sin  0dcf>  +  e  dt  (3.5) 

where  A  and  v  are  functions  of  r  only.  The  corresponding  energy-momentum 
tensor  is  given  by 


where  p  is  the  density  of  the  cosmological  fluid  and  p  is  the  pressure. 
Calculating  Christoff el  symbols,  Ricci  tensor  and  scalar  curvature  in 
the  usual  way  one  can  then  compute  the  field  equations  (3.4).  Tolman  [12] 
shows  that  one  of  these  equations  reduces  to 


* 
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(p  +  p)v'  =  0  ,  (3.7) 

which  implies  that  either  (i)  v'  =  0  or  (ii)  p  +  p  =  0  or  (iii)  both 
v'  and  p  4-  p  are  zero.  Case  (i)  yields  the  Einstein  static  universe, 
case  (ii)  gives  the  de  Sitter  static,  universe  and  case  (iii)  gives  the 
special  relativity  line  element.  These  three  cases  exhaust  all 
possibilities  for  a  spherically  symmetric  static  universe.  Let  us 
investigate  the  cases  in  some  detail. 

Case  (i) :  v'  =  0  .  Simple  mathematical  manipulation  yields  the  line 
element 

,  2  dr2  2  , „ 2  2  .  2  _,2  .  ,.2 

ds  =  -  — '  -  -  -  -  •  -  r  d6  -  r  sm  0d<}>  +  dt  (3.8) 

1-r  /R 

which  is  a  standard  form  of  the  metric  for  the  Einstein  model.  The 
constant  R  in  (3.8)  is  defined  by 

A  -  8irp  =  ~  .  (3.9) 

R 


The  coordinate  substitution  r  =  R  sin  ip  yields  the  metric  form 

ds^  =  -  R^  (d\p^  +  sin^  xpd0 ^  +  sin^  ip  sin^  6dcp ^ )  +  dt^  .  (3.10) 


Introducing  an  extra  variable  by  means  of  the  equations 


2  2  1/2 

y  =  R(1  -  r  /R  )  ,  y  =  r  sin  0  cos  cj>  , 


y^  =  r  sin  0  sin  cp  ,  y^  =  r  cos  0  , 


}  (3.11) 


the  line  element  can  be  written  in  the  form 


2  2  2  2  2  2 
ds  =  -  dy1  -  dy2  -  dy3  -  dy4  +  dt 


(3.12) 
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However  equations  (3.11)  imply  that 


yl  +  y2  +  y3  +  y4  =  R 


(3.13) 


Equation  (3.12)  leads  to  the  conclusion  that  our  original  four-dimensional 
space- time  is  embedded  in  a  five-dimensional  Euclidean  space.  From 
equation  (3.13)  we  conclude  that  the  three-dimensional  hypersurface  which 
composes  the  spatial  region  of  the  Einstein  universe  is  in  fact  the 
surface  of  a  four-dimensional  hypersphere  with  radius  R.  This  justifies 
calling  R  the  radius  of  the  universe. 


In  accordance  with  the  metric  form  (3.10)  the  spatial  volume 

of  the  universe  is  given  by 

2tt  ti  ir 


v 


”  0  J  0 


3  2  2  3 

R  sin  ip  sin  0dij;d9d<f)  =  2tt  R 

J0 


(3.14) 


and  the  circumference  is  2ttR  =  An  interesting  feature  of  the  Einstein 

model  is  that  the  surface  area  of  a  sphere  of  radius  r  is  less  than 
2 

4?rr  .  Einstein  [13]  gives  an  excellent  qualitative  discussion  of  the 
geometry  of  this  universe. 

A  disappointing  feature  of  the  model  is  the  absence  of  an 
intrinsic  red-shift  in  the  light  from  the  nebulae.  Investigation  of  the 
geodesic  equations  corresponding  to  the  metric  (3.8)  yields  the  result 
that  a  free  particle  which  is  at  rest  with  respect  to  the  spatial 
coordinate  at  any  time  must  in  fact  be  permanently  at  rest.  Hence  there 
can  be  no  Doppler  effect  on  the  light  received  by  an  observer  at  the 
origin  r  =  0  from  a  nebula  at  r  =  r^  where  both  observer  and  light 
source  are  assumed  to  be  at  rest  relative  to  the  matter  in  their  respective 
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neighbourhoods . 


Case  (ii) :  p  +  p  =  0  .  Physical  arguments  lead  us  to  conclude  that 
p  =  p  =  0  i.e.  space  is  completely  devoid  of  matter  or  radiation. 
Evidently  such  a  model  is  unsuitable  as  a  picture  of  the  actual  universe; 
however  it  does  have  an  intrinsic  red-shift  in  the  light  received  from 
test  particles  of  negligible  mass  and  so  it  has  been  of  great  interest 
to  cosmologists .  The  line  element  for  this  de  Sitter  universe  can  be 
derived  in  the  form 

,2  dr2  2  _2  2  .  2  „  ,  ,2  ,  r2N  ,  2  _  „  _s 

ds  =  -  — —  -  r  d0  -  r  sin  0d4>  +  (1  -  — )  dt  (3.15) 

1-r  /R  RZ 

where  ■—  5  ^  „  If  R  is  regarded  as  the  radius  of  the  universe  then 

the  model  is  either  spatially  closed,  open  and  flat  (special  relativity) 
or  open  and  curved  according  as  A  is  positive,  zero  or  negative. 

Evidently  metric  (3.15)  can  be  obtained  by  setting  m  =  0  in 
Schwarzschild' s  exterior  solution  for  the  gravitational  field  surrounding 
a  spherical  body  of  mass  m. 


Since  a  light  ray  travels  on  a  null  geodesic  the  velocity  of 
light  travelling  from  or  towards  the  origin  is  given  by 


dr 

dt 


(3.16) 


Integrating  (3.16)  between  the  limits  r  =  0  and  r  =  R  we  find  that 
it  would  take  an  infinite  time  interval,  as  measured  by  an  observer  at 
the  origin,  for  a  light  ray  to  travel  from  the  origin  to  r  =  R  . 

Hence  this  distance  is  called  the  observer's  "horizon  of  the  universe" 
as  he  can  never  obtain  information  concerning  events  occurring  at  or 


. 


22 


beyond  it.  The  geodesic  equations  yield  the  orbits  of  test  particles 
and  combining  them  with  (3.16)  one  can  compute  as  an  expression  for 
red-shift  in  the  waveleng-th  of  light  received  from  a  particle 


A  +  AX 

A 


k 


2  2 

1  -  r  /R 


+ 


k 


2  2  2 
(1  -  r  /R  ) 


dr 

dt 


(3.17) 


where  k  is  a  positive  constant  and  is  the  radial  velocity  of  the 

particle  at  the  time  of  emission.  Thus  we  see  that  in  the  de  Sitter 
universe  one  would  expect  to  observe  a  wavelength  shift  in  the  light 
from  freely  moving  particles.  Coupling  this  result  with  Weyl’s  hypothesis 
it  is  possible  to  show  that  the  relation  between  red-shift  and  distance 
is  approximately  linear,  as  is  suggested  by  observational  data. 


This  account  of  the  Einstein  and  de  Sitter  static  models  leads 
us  to  conclude  that  neither  can  be  regarded  as  a  good  description  of  the 
actual  universe.  Hence  it  is  necessary  to  consider  a  less  restrictive 
line  element  in  constructing  new  models. 

In  chapter  11  we  considered  the  cosmological  principle  which 
leads  to  the  Robertson-Walker  metric.  The  first  assumption  was  that 
there  exists  a  cosmic  time  t  which  is  orthogonal  to  the  three- 
dimensional  hypersurfaces  describing  the  spatial  properties  of  the 
universe.  This  condition  enables  us  to  construct  a  system  of  co-moving 
coordinates  in  which,  at  any  moment,  the  matter  of  the  universe  is  on 
the  average  at  rest.  Further  the  line  element  can  be  written  in  the  form 

,2  ,  i  ,  j  ,  ,  2 

ds  =  -  g . .  dx  dx  +  dt 

ij 


(3.18) 


' 
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where  i,j  =  1,2,3  and  the  g_  are  functions  of  x  ,x  ,x  ,t.  Hence  this 

line  element  is  non-static.  However  the  assumption  of  isotropy  implies 
that  the  spatial  coordinates  again  appear  in  spherically  symmetric  form. 
Thus  (3.8)  can  be  written 

,  2  G(r,t)  2  2  ,Qs 

ds  =  -  e  ’  dl  +  dt  (3.19) 


where 

2  2  2  2  2  2  2 
dl  =  dr  +  r  d0  +  r  sin  0dcf) 


(3.20) 


and  G(r,t)  is  independent  of  0  and  (f)  as  a  result  of  isotropy.  The  third 
assumption  concerning  the  equivalence  of  any  two  points  on  a  hypersurface 
t  =  constant  can  be  shown  to  lead  to  the  result 

G (r , t)  =  g(t)  +  f (r)  ,  (3.21) 


where  g  and  f  are  arbitrary  functions  of  t  and  r  respectively.  The 
metric  (3.19)  now  reduces  to  the  form 


ds 


2 


,g(t)  + 


f(r)di2 


+  dt‘ 


(3.22) 


Investigation  of  the  cosmological  field  equations  for  this  line  element 
enables  us  to  determine  f (r)  explicitly  and  we  finally  obtain  the 
Robertson-Walker  line  element 


ds2  = 


g(t) 


[1  +  !i  £-  ]2 

L  4  2J 

ro 


9  2  2  2  2  2  2 

(dr  +  r  d0  +  r  sin  0d<f>  )  +  dt 


(3.23) 


where  k  =  0,  +1,  -  1  according  as  the  curvature  of  space-time  is  zero, 
positive  or  negative,  and  r^  is  a  constant.  See  Adler,  Bazin  and 


' 
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Schiffer  [14]  for  a  detailed  treatment  of  this  derivation. 


In  astronomy  it  is  often  convenient  to  define  a  new  radial 
coordinate  u  =  —  and  a  new  function 

ro 


R2 (t)  =  rQ2  e8(t) 


(3.24) 


Substituting  in  (3.23)  we  arrive  at  the  form 


ds‘ 


2  2  2  2  2  2 
„2,.x  du  +  u  d0  +  u  sin  0dd>  .  „  2 

n  +  k  2,2  +  dt 

(X  +  4  u  ) 


(3.25) 


R(t)  is  called  the  radius  of  the  universe.  As  in  the  Einstein  model  we 
can  show  that  the  space  part  of  (3.25)  is  the  line  element  of  a  three- 
dimensional  hyperspherical  surface  embedded  in  a  Euclidean  4-space.  There 
are  three  cases  to  consider: 

(i)  k  =  1,  then  R  is  real  and  the  hypersphere  has  uniform  positive 

curvature.  The  space  is  closed  and  the  surface  area  of  a  sphere  of 

2 

radius  r  is  less  than  4TTr 


(ii)  k  =  0,  then  R  is  infinite  i.e.  curvature  is  zero  and  the  3-space 
is  Euclidean. 


(iii)  k  =  -1,  then  R  is  imaginary  and  the  space  is  a  "pseudohypersphere. 11 

It  is  an  open  space  with  uniform  negative  curvature  in  which  the  surface 

2 

area  of  a  sphere  of  radius  r  is  greater  than  4iTr 

Finally,  it  is  not  difficult  to  demonstrate  that  the  expression 
for  wavelength  shift  in  a  Robertson-Walker  universe  is  given  by 


«  -  IJ  -o  I*  -  ,0  J.  Sfi 
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(3.26) 


which  is  identical  with  the  Newtonian  result  (2.31)  and  where  the  terms 
have  the  same  interpretation.  Once  again  it  is  evident  that  the  condition 
for  the  shift  to  be  towards  the  red  is  that  >  R^  i.e.  the  universe  is 

expanding.  Doppler's  law,  =  v  ,  implies  that 

A 


z  = 


AX 

X 


^2 


(3.27) 


where  L  is  the  proper  distance  between  observer  and  light  source. 
Comparing  (3.27)  with  Hubble's  Law 


AX 

X 


L 

T 


(3.28) 


where  T  is  Hubble's  constant,  we  get  as  a  first  approximation  for  T 


1 

T 


(3.29) 


Having  derived  the  Robertson-Walker  metric  and  demonstrated  that  it 
explains  red-shift  we  are  now  ready  to  discuss  the  various  non-static 
model  universes  to  which  it  gives  rise. 
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CHAPTER  IV 


Non-static  Universes 


The  field  equations  (3.4)  for  the  metric  (3.25)  are  given  by 


8ttp  =  — 


3k  3R* 2 

2  +  2 
R  R 


-  A 


(4.1) 


8np  = 


k_  2R"'1 

R2  R2  R 


+  A 


(4.2) 


where  we  have  chosen  a  co-moving  coordinate  system  such  that 


dr  _  d_0_  _  _d<j>_  _  „  _dt  _ 

ds  ds  ds  *  ds 


(4.3) 


and  the  energy-momentum  tensor  is  given  by  (3.6).  Elementary  algebraic 
manipulation  of  (4.1)  and  (4.2)  leads  to  the  equivalent  system 


4tt  (p  +  3p)  =  A  - 


3R" 

R 


(4.4) 


4tt  (p  +  p)  =  —  + 

R 


k  .  R’2  -  RR"  k  d  rR’ 


R 


k_  _  a_  ,-k.  a 
2  dt  R  J 


R 


(4.5) 


Differentiating  equation  (4.1)  with  respect  to  time  and  combining  the 
result  with  (4.5)  we  get 

R  +  3(p  +  p)R’  =  0  .  (4.6) 

2 

Multiplying  (4.6)  by  the  integrating  factor  R  and  rearranging  we  get 

d  ,  _3.  d  ,  3. 

—  (PR  }  +  P  (R  )  - 


0 


(4.7) 
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In  accordance  with  the  Robertson-Walker  line  element  (3.25)  it  is  clear 

3  3 

that  proper  volume  is  proportional  to  R  and  consequently  pR  is 
proportional  to  the  corresponding  mass.  Hence  if  V(t)  is  the  volume 
and  M  =  pV  is  the  mass  equation  (4.7)  can  be  written 

dM  +  pdV  =  0  (4.8) 

Since  M  measures  the  energy  content  of  the  volume  element  V  and  pdV 
is  a  measure  of  the  work  done  against  pressure  (4.8)  is  the  law  of 
conservation  of  energy  (mass) .  Note  that  the  Newtonian  counterpart 
of  (4.8)  is  simply  dM  =  0  ,  i.e.  mass  alone  is  strictly  conserved. 

In  the  special  case  R  constant,  equations  (4.1)  and  (4.2) 

become 

8rp  =  -  A  (4.9) 

R 

8ttP  =  A  -  —  (4.10) 

R 

and  if  k  =  1  these  are  the  equations  for  pressure  and  density  in  the 
Einstein  static  universe.  Assuming  p  and  p  are  positive  quantities 
equations  (4.9)  and  (4.10)  imply  that 

—  <  A  <  .  (4.11) 

R  R 

Hence  it  follows  that  A  has  the  order  of  magnitude  of  the  inverse 
square  of  the  radius  of  the  universe,  which  explains  why  it  cannot  be 
expected  to  have  any  significance  on  the  scale  of  the  solar  system. 
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In  constructing  non-static  models  we  can  arbitrarily  specify  the 
values  of  A  and  k,  and  the  equation  of  state  of  the  fluid,  p  =  p(p)  . 

Once  the  equation  of  state  is  prescribed  (4.7)  can  be  written  in  the  form 


dR  1 r  dp 
R  3  '•p  +  p(p) 


In  the  simple  case  p  =  0  (4.12)  yields 


R  p  =  constant 


(4.12) 


(4.13) 


If  the  pressure  is  due  entirely  to  radiation  the  equation  of  state  is 
p  =  ^p  and  integration  of  (4.12)  implies 

4 

R  p  =  constant  .  (4.14) 

Friedmann  in  1922  found  a  non-static  solution  with  A  =  0  . 

This  was  the  first  physically  reasonable  cosmological  solution  of  the 
field  equations  without  A  and  its  discovery  led  Einstein  to  discard  the 
cosmological  term  which  he  felt  to  be  aesthetically  undesirable.  The 
original  Friedmann  solution  corresponds  to  the  case  A  =  k  =  0.  Prompted 
by  the  physical  evidence  that  the  pressure  is  very  much  smaller  in 
magnitude  than  the  density  in  the  actual  universe,  the  simplification 
p  =  0  is  also  made.  Equation  (4.2)  then  becomes 

R'2  +  2RR"  =  0  ,  (4.15) 

which  is  easily  integrated  to  give 

2 

R’  R  =  C  ,  (C  arbitrary  constant).  (4.16) 


Integrating  (4.16)  between  appropriate  limits  we  obtain 


‘ 
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3/2 


■r> 


(4.17) 


It  is  always  possible  to  choose  coordinates  such  that  R  =  0  when 
t  =  0  so  that  (4.17)  can  be  written  in  the  simpler  form 


R 


(2)3/2  cl/3  t 2/3 


(4.18) 


This  universe  has  already  been  derived  in  our  discussion  of  Newtonian 
cosmology.  It  is  a  monotonically  expanding  universe  which  starts  to 
expand  from  a  singular  state  at  time  t  =  0.  Such  a  model  is  sometimes 
called  "exploding".  Equation  (4.1)  becomes 

3Rl  ^ 

8ttP  =  — r-  (4.19) 

R 

and  comparing  this  with  equation  (3.29)  we  get 

8ttP  =  -2-  (4.20) 

T 


Taking  Hubble's  constant  T  to  have  the  value  T  =  (4.1  1  2)10^  sec.  we 
obtain  for  the  present  density  of  the  universe 


.  _29  3 

p  *  (1.1  +  0.5)10  g/cm 


(4.21) 


-27  -31  3 

which  lies  within  the  experimental  range  of  10  to  10  g/cm 


At  the  present  time,  with  no  observational  evidence  to  the 
contrary,  A  =  0  seems  to  be  the  most  reasonable  hypothesis.  Friedmann 
investigated  more  general  solutions  by  putting  A  =  0  but  imposing  no 
restriction  on  k.  The  field  equations  (4.1)  and  (4.4)  become, 
respectively, 


. 
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8-rrp 


(4.22) 


and 


% 

4tt(p  4-  3p) 


3R" 

R 


(4.23) 


If  we  accept  the  physically  reasonable  assertion  that  p  and  p  must  be 
non-negative,  equation  (4.21)  implies  that  R"  <  0  ,  hence  the  rate  of 
expansion  must  be  retarded  and  from  elementary  calculus  we  conclude  that 
R  can  never  have  a  minimum  or  a  point  of  inflection.  For  simplicity  we 
again  assume  that  pressure  is  negligible;  then  by  eliminating  p  between 
equations  (4.22)  and  (4.23)  with  p  =  0  ,  R  can  be  expressed  as  a  function 
of  t  by  the  differential  equation 


2RR"  +  R'2  +  k  =  0 

Rewriting  (4.24)  in  the  equivalent  form 

^(RR'2)  +  kf  -  0 

integration  immediately  yields 

2  R  -  kR 

R*  =  S  -  (R  constant) 

R  s 

Inserting  (4.26)  into  (4.24)  we  obtain 

R 

....  s 


(4.24) 


(4.25) 


(4.26) 


(4.27) 


Now  this  result  must  be  compatible  with  the  requirement  that  R"  is 

negative,  hence  R  is  a  positive  constant.  Equations  (4.23)  and  (4.27) 
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with  p  =  0,  imply  that 

Rg  =  2  .  (4.28) 

However  the  quantity  in  parentheses  represents  the  total  mass  M  of  the 
universe  so  that  (4.28)  becomes 

R  =  2M  .  (4.29) 

s 


This  is  the  radius  of  discontinuity  of  the  Schwarzschild  exterior  solution. 
We  must  consider  the  three  cases  k  =  0,  +  1,  -  1. 


Case  (i) :  k  =  1.  Then  equation  (4.26)  becomes 


RR’ 


2 


R  -  R 
s 


(4.30) 


and  the  substitution  R  =  R  sin  u(t)  permits  us  to  solve  (4.30)  in  the 

s 

parametric  form 


t  =  ^R  (2u  -  sin  2u)  (4.31) 

2  s 

R  =  t;R  (1  -  cos  2u)  .  (4.32) 

2  s 


Equations  (4.31),  (4.32)  are  recognised  as  the  equations  of  a  cycloid. 

When  t  =  u  =  0  we  get  R  =  0;  then  as  t  increases  the  radius  R  expands  up 

to  a  maximum  value  R  at  time  t  =  tS-ttR  after  which  R  contracts.  A 

s  2  s 

universe  of  this  type  is  called  "oscillating." 


Case  (ii) :  k  =  0.  It  has  already  been  shown  that  A  =  k  =  0  leads  to 

the  monotonically  expanding  universe  described  by  equation  (4.18). 


. 
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Case  (iii) :  k  =  -  1.  Then  equation  (4.26)  becomes 

2  R 

R'  =  “jf  +  1  (4.33) 

2 

and  the  substitution  R  =  R  sinh  u(t)  yields  the  solution 

s 

t  =  ^R  (sinh  2u  -  2u)  (4.34) 

l  s 

R  =  ^R  (cosh  2u  -  1)  .  (4.35) 

l  s 

When  u  =  t  =  0  we  get  R  =  0  after  which  R  increases  steadily  to  infinity. 
Returning  to  equation  (4.33)  we  see  that  for  large  R 

R'  -  1  ,  R  -  t  (4.36) 

and  recalling  that  the  velocity  of  light  is  unity  it  follows  that  the 
velocity  of  expansion  is  approximately  the  velocity  of  light.  It  is 
also  interesting  to  observe  that  R  =  0  implies  R'  =  00  and  so  the  universe 
has  an  explosive  birth. 

The  behaviour  of  these  models  appears  to  be  physically 
reasonable  except  for  the  singular  state  which  exists  at  the  moment  of 
"birth."  The  singularity  is  a  result  of  neglecting  the  pressure  term 
in  the  field  equations  because  when  the  radius  of  the  universe  is  small 
and  its  matter  and  radiation  are  in  a  highly  condensed  state  it  is 
reasonable  to  expect  a  pressure  term  of  large  magnitude. 

Theoretically  it  should  be  possible  when  p  can  be  experimentally 
established  to  within  a  sufficient  degree  of  accuracy  to  determine  which 
of  the  above  three-eases- corresponds  to  the- actual  universe.  However 
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with  the  existing  state  of  observational,  data  one  can  only  speculate  as 
to  which  type  of  universe  is  aesthetically  more  desirable  since  all  are 
mathematically  valid. 

If  for  the  sake  of  permitting  a  greater  multiplicity  of  models 
we  retain  A  in  the  field  equations  then  it  is  necessary  to  consider 
equations  (4.1)  and  (4.2)  in  the  general  case.  With  the  present  state 
of  experimental  cosmology  there  is  no  evidence  to  indicate  a  preference 
for  including  the  cosmological  term,  though  it  is  of  course  satisfying 
to  investigate  all  the  possibilities  for  a  mathematical  world  model. 

A  particular  solution  with  A  may  be  obtained  by  assuming  that  the  rate 
of  expansion  of  the  universe  is  given  exactly  by  equation  (3.29),  viz. 

t"  “  i  (4.37) 


where  T  is  Hubble’s  constant.  Integration  yields 


R  Roe 


t/T 


Rq  constant. 


(4.38) 


Substituting  in  (4.4)  and  (4.5)  we  obtain  the  result 

p  =  p  -  k  =  0  ,  (4.39) 

i.e.  the  universe  is  completely  empty.  However  it  supports  an  intrinsic 
red-shift  in  the  light  received  from  the  nebulae  (particles  of  negligible 
mass) .  But  these  are  precisely  the  properties  of  the  de  Sitter  static 
model  described  in  chapter  111.  Indeed  Tolman  [12]  defines  a  coordinate 
transformation  which  exhibits  the  equivalence  of  the  respective  metrics. 
This  universe  has  the  desirable  property  that  it  has  no  singularity 
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i.e.  R  is  neither  zero  nor  infinite  for  a  finite  value  of  t. 
Furthermore  it  does  not  have  a  "birth"  or  "death"  ,  a  consequence  of 
its  infinite  time-scale  and  lack  of  singularity.  It  is  interesting 
to  note  that  all  expanding  universes  with  a  constant  amount  of  matter 
and  radiation  must  ultimately  degenerate  into  a  de  Sitter  universe. 

It  is  sometimes  desirable  to  specify  a  certain  property  in 
advance  and  then  to  investigate  the  nature  of  the  universe  which 
results  from  the  given  assumption.  Lemaitre  [15],  for  instance,  in 
1927  examined  the  case  of  closed  models  (k  =  1)  containing  incoherent 
matter  (p  =  0)  and  radiation  (p  =  ^p  ,  p  is  energy  density) . 

The  energy  equation  (4.7)  becomes 

{(pm  +  3p)R3}  +  p  -f^(R3)  =  0  (4.40) 


where  pm  is  matter  density  and  the  total  density  is  p 
the  law  of  conservation  of  matter  implies  that 

3  4 

PmR  =  constant,  pR  =  constant 

from  which  we  conclude  that 

8ttp  =  — 

“  R3 

8ttP  =  ^ 

R 

8irp  =  8tt  (p  +  3p)  =  +  ^f 

R  R 


m 


+  p 
Ke 


Now 


(4.41) 


}  (4.42) 


J 


where  a  and  8  are  constants. 


Simple  manipulation  of  (4.1)  and  (4.42) 


now  yields 


' 
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dR 

dt 


±  t 


3  a/2 


(4.43) 


3 


R 


Integration  of  (4.43)  has  been  discussed  by  de  Sitter  [16].  Note  that 
in  the  special  case  A  =  6  =  0  equation  (4.43)  becomes  equation  (4.30) 


where  R  =  ^ 
s  3 


In  this  chapter  our  discussion  has  been  primarily  concerned 


with  models  for  which  A  =  0  .  Inclusion  of  A  in  the  field  equations 
leads  to  many  other  models,  the  de  Sitter  empty  universe  being  a 
particularly  interesting  example.  However  integration  of  the  formulae 
for  R  in  terms  of  t  is  greatly  complicated  by  the  presence  of  terms 
involving  A  and  as  has  been  stated  above  there  is  at  the  present  time 
no  physical  argument  in  favour  of  retaining  the  cosmological  term  in 
the  field  equations. 
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CHAPTER  V 

The  Godel  Solution :  A  Rotating  Universe . 

Godel’s  universe  is  of  historical  interest  in  that  it  is  the 
earliest  relativistic  model  which  does  not  have  a  Robertson-Walker  metric. 
In  1949  Godel  [11]  demonstrated  that  the  line  element 


dsZ  =  -  (dx')2  +  \  e2aX’(dx2)2  -  (dx3)2  +  2eaX'dx2dx4  +  (dx4)2 


(5.1) 


where  a  is  a  constant  yields  a  solution  of  the  field  equations  when  the 
cosmological  fluid  is  idealized  as  an  incoherent  matter  distribution.  The 


metric  tensor  g . .  and  its  inverse  g1^  are  given  by 
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(5.2) 


Computing  the  Ricci  tensor  in  the  usual  manner  we  obtain 
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and 


R 


•% 


(5.4) 


The  energy-momentum  tensor  for  an  incoherent  matter  distribution 
at  rest  in  a  co-moving  coordinate  system  is  given  by 


,ij  dx1  dx^  -i  j 

-  p  dT  dT  “  p64  S4 


(5.5) 


where 


&1.  is  the  Kronecker  delta  from  which  it  follows  that 
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T.  . 
1J 
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(5.6) 


Comparing  (5.3)  and  (5.6)  we  get 

2 

R.  .  =  -  —  T.  . 

ij  P  iJ 


(5.7) 


Substituting  (5.4)  and  (5.7)  into  Einstein's  field  equations  (3.4)  we 
obtain 


+  8tt]  T 


ij 


+ 


(5.8) 


which  is  identically  satisfied  if 


1  2 
2 


a 


(5.9) 


Note  that  if  a  =  0  we  get  A  =  p  =  0  which  implies  that  space  is  flat. 
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Thus  a  may  be  regarded  as  a  parameter  which  measures  the  curvature  of  space, 

Let  us  now  recall  that  an  Einstein  static  universe  with  an 
incoherent  matter  distribution  has  the  energy-momentum  tensor  (5.5). 

This  leads  to  the  surprising  conclusion  that  the  same  T  gives  rise 
to  two  fundamentally  different  universes.  This  result  substantiates 
the  argument  that  Mach’s  principle  is  not  an  intrinsic  property  of 
general  relativity  since  one  can  no  longer  be  sure  that  a  given  T 
uniquely  specifies  an  appropriate  metric. 


An  interesting  property  of  the  Godel  solution  is  that  it 
does  not  have  a  cosmic  time,  as  do  earlier  Robertson-Walker  universes. 
It  can  be  shown  that  the  world  lines  of  matter  at  rest  in  the 
co-moving  coordinate  system  cannot  be  everywhere  orthogonal  to  a  one- 
parameter  family  F  of  three-dimensional  hypersurf aces .  A  covariant 
necessary  condition  that  a  vector  field  vr  be  everywhere  orthogonal  to 
F  is 


Vr  V  i  =  0 

[r  s,t] 


(5.10) 


where  the  tensor  vr  v  .  ,  is  the  skew-symmetric  part  of  v  v  and 

lr  s,tj  rs,t 

the  comma  denotes  partial  differentiation. 

The  vector  which  specifies  matter  at  rest  in  the  Godel  model 
is 


v1  =  (0,  0,  0,  1)  ,  vr  =  (0,  e“X  ,  0,  1)  ;  (5.11) 
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for  this  vector 


v 

r,s 


1 

ax 

ae 


for  r  =  2,  s  =  1 
otherwise 


Hence  the  tensor  vr  v  ,  is  given  by 

[r  s,t] 


(5.12) 
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E  ±  6ijk 

3!  rst 

n  1 

I  1  ax 


even  permutation  of  4 

odd  permutation  of  4, 
otherwise 


1,  2 

1,  2 


(5.13) 


i.e.  vr  v  i  is  not  identically  zero  which  implies  that  there  is  no 
[r  s,t] 

vector  field  which  is  everywhere  orthogonal  to  F. 


In  order  to  gain  an  intuitive  understanding  of  the  rotating 
nature  of  this  model  we  note  that  the  line  element  describing  a  flat 
space  with  cylindrical  coordinates  r,<p  and  z  rotating  about  the  axis  of 
z  with  constant  angular  velocity  w  can  be  written. 

2  2222  2  222 
ds  =  -  dr  -  r  d<f>  -  dz  -  2cor  dcf>d t  +  (l-oo  r  )dt  .  (5.14) 

Comparing  this  with  Godel's  metric  (5.1)  we  observe  the  similarity  of 

the  two  forms  and  note  in  particular  that  each  has  a  cross  term 

involving  the  time  coordinate.  Hence  we  may  tentatively  identify  x^ 

2  3 

with  r,  x  with  <j>  and  x  with  z. 

Now  consider  the  vector  ft1  defined  by 
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1  e 

n  = 


xrst 


-  v  v 

/T~  [r  s,t] 


(5.15) 


irst 


where  is  the  well-known  permutation  symbol.  Computing  the  components 

of  ^  in  a  flat  space  for  which  /-g  =  1  the  vector  (£t  ,  associated 

/  .  v  4 

with  (— v  ,  v  )  is  seen  to  be 


^  =  +  v^(V  x  v)  +  (v  x  +  (v  x  v)v^ 


=  +  v* (V  x  v) 


(5.16) 


Let  us  now  associate  v  with  the  vector  field  describing  uniform  counter¬ 
clockwise  rotation  about  the  z-axis  of  a  Cartesian  3-space.  Suppose 
co  is  the  angular  velocity  of  the  system.  Clearly 


v  =  (coy,  -cox,  0) 


v4  =  1  ’ 


(5.17) 


from  which  it  follows  that 


=  (0,  0,-2co,  0) 


(5.18) 


For  the  Godel  metric  can  be  computed  directly  from  (5.13)  and  (5.15) 


to  give 


ft1  =  (0,  0,-/2a,  0) 


(5.19) 


Hence  by  analogy  with  (5.18)  we  may  conclude  that  (5.19)  implies  that 
the  co-moving  matter  of  the  Godel  universe  possesses  a  constant  angular 
velocity 


a 

'  h 


co 


5 


(5.20) 


. 
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about  the  axis  of  x~ .  It  is  interesting  to  note  that  we  now  have  a 
matter  distribution  at  rest  relative  to  a  co-moving  coordinate  system 
which  is  non- inertial .  This  is  clearly  inconsistent  with  Mach’s 
principle  as  stated  above. 

The  motion  of  a  test  particle  can  be  considered  by  investigating 
the  geodesic  equations  for  the  Godel  metric.  One  arrives  at  the  conclusion 
that  a  test  particle  projected  radially  from  the  origin  of  coordinates 
toward  a  distant  nebula  will  not  travel  in  a  straight  line  towards  its 
goal,  as  Mach's  principle  would  lead  us  to  expect.  Instead  the  particle 
will  spiral  outward;  it  is  relative  to  the  path  of  such  a  particle  that 
the  matter  of  the  universe  rotates.  The  tangent  of  this  path  is  called 
the  "compass  of  inertia".  Hence  we  may  say  that  the  matter  of  a  Godel 
universe  rotates  relative  to  the  compass  of  inertia,  or  vice  versa. 

More  general  line  elements  describing  rotating  universes  are 
considered  by  Synge  [17]  and  Wright  [18],  The  most  unsatisfactory 
feature  of  the  Godel  solution  is  the  absence  of  a  red-shift  which  renders 
it  unsuitable  as  a  model  of  the  actual  universe.  In  order  to  obtain  a 
Godel  type  universe  which  is  compatible  with  observational  data  it  is 
necessary  to  combine  expansion  with  rotation. 


. 

. 
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CHAPTER  VI 

Some  New  Solutions  of  the  Cosmological  Field  Equations. 


form 


In  this  chapter  we  shall  investigate  the  field  equations  in  the 


R.  . 
1J 


2Sij R  +  Agij 


T.  . 
13 


(6.1) 


where  the  metric  tensor  is  assumed  to  be  given  by 


/2a  1  0  ° 

10  0  0 
0  0  e“2T  0 


VO  0  0  e  2  V 


(6.2) 


A  more  general  metric  has  been  considered  by  Wyman  and  Trollope  [19],  [20] 
in  connection  with  electromagnetic  null  fields  in  general  relativity  and 
their  success  in  solving  the  Einstein-Maxwell  equations  suggests  that  a 
similar  metric  may  also  yield  some  interesting  cosmological  solutions. 
However  their  original  g_  leads  to  rather  unwieldy  field  equations  and  so 
it  has  been  modified;  consequently  the  present  analysis  is  restricted  to 

the  form  (6.2),  in  which  a  and  x  are  functions  of  the  four  coordinates 

1  2  3  4  3  4 

x  ,  x  ,  x  ,  x  .  The  coordinates  x  and  x  are  assumed  to  be  pure  space 

1  2 

coordinates  whereas  x  and  x  are  null  coordinates.  In  order  to  separate 
the  space  and  time  elements  consider  the  transformation 


-2  -1  -2  -1 

lx-x  2  x  +  x  3-3  4-4 

x  =  -  ,  x  =  -  ,  x  =  x  ,  x  =  x 

/2  /2 


(6.3) 


•  • 


'TTO  -3  '  '■  is  79i;  2003  8309109X3  9Clt3  bfTB  -ODBC  I  »fl3 
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In  terms  of  the  new  coordinates  we  have 


8ij 


N 

P 

1 

H* 

-a, 

0 

0  ^ 

-a. 

a+1 

0 

0 

0 

0 

-2t 

e 

0  | 

V  o 

0 

0 

e-2r) 

(6.4) 


If  it  is  assumed  that  a  is  less  than  one  then  is  a  timelike  coordinate. 

If  we  further  require  that  x  and  a  vanish  at  the  origin  then  g  reduces 
to  the  special  relativity  metric  with  signature  +2.  Although  the  field 
equations  will  be  calculated  with  respect  to  (6.2)  nonetheless  it  will 
frequently  prove  convenient  to  refer  to  the  barred  coordinate  system  in 
order  to  gain  physical  insight.  The  presence  of  the  off-diagonal  terms 
in  (6,4)  leads  to  the  hope  that  the  solutions  may  correspond  to  a  rotating 
universe  of  a  nature  similar  to  that  of  Godel. 


Let  us  assume  that  the  universe  is  filled  with  dust  for  which 
we  would  have  p  =  0.  The  energy-momentum  tensor  is  then  defined  by 


TlJ  =  puV 


(6.5) 


i  dx 

where  p  is  proper  density  and  u  =  — -  is  the  four-velocity  characterizing 

Cl  s 

the  flow  of  the  cosmological  fluid.  Suppose  now  that  the  barred  coordinate 
system  is  co-moving,  then 


u1  =  TT-  =  (  —  ,0,0,0) 

ds 


(6.6) 


It  follows  from  the  transformation  law  for  a  contravariant  vector  that 


•  ,/ 


. 
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i  9x"  -i  1  ,  „ 

u  =  ~tt  u  =  ~nrm  °»  o) 

3xJ 


/2 (1-a) 


(6.7) 


The  covariant  velocity  vector  is  given  by 


— (-2a  +  1,-1,  0,  0) 
/2 (1-a) 


(6.8) 


and  the  covariant  energy-momentum  tensor  T..  =  pu.u. 

ij  i  J 


obtained  in  the  form 


T.  .  =  p 
ij 


(  (l-2a) 2 
2 (1-a) 


2a-l 


2 (1-a) 
0 


0 


2a-l 
2 (1-a) 

1 

2 (1-a) 
0 
0 


0  0 


0  0 


0 

0 


\ 


0 

0 


J 


(6.9) 


The  scalar  T  =  g^T..  can  be  shown  to  be 

ij 


T  =  -  P 


(6.10) 


and  if  we  now  multiply  the  field  equation  (6.1)  throughout  by  g1^  we 
obtain 


R  =  4A  -  T 


(6.11) 


In  view  of  equation  (6.10)  we  get 


R  =  4A  +  p 


(6.12) 


The  field  equations  can  now  be  written  in  the  convenient  form 


R.  . 
ij 


i  g  (p  +  2A) 


(6.13) 
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The  formula 


{jk}  2  g  (8rj,k+8rk,j  8jk,r) 


yields  the  following  non-zero  Christoff el  symbols: 


{ii}  ■ 


r  1X  r  1 \  _  ~2X 

“  01  „  »  1  -J  O  !  ^  T  T e  » 


,2 


33  L44 


{11}  -  “,1  +  2aa,2  ’  {12}  =  “,2 


f13}  “,3  ’  (14}  “,4 


*33*  t44t  <\l  2aT,2)  e 


<11}  = 


-2t 

aj3e 


{  3}  =  {  ^}  = 
^13  ^  l14^ 


-  T 


,1  ’ 


{  3}  =  {  4}  = 
X23J  z4; 


-  T 


,2  ’ 


{  3)  =  {  = 
03^  ^34; 


-  T 


,3  ’ 


{34> 


,4  * 


{  3> 

^44 ; 


=  T , 3  ’  {11}  =  '  a,4e 


-2x 


{33}  T , 4 


(6.14) 


>  (6.15) 


) 
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Computing  the  Ricci  tensor  R. .  in  the  usual  manner  and  substituting 

J-j 

in  (6.13)  we  obtain  the  following  system  of  differential  equations: 


-  2x  ^  -  2ota  22  +  +  2(x  -  2a  0t  ^  +  2a  ^ 


+  4aa  0x  =  2aA  -  ap  +  ,  , 

,2  ,2  2(l-a) 


"a,22  2ts12  4  2t,2  (cts2  4  T,r 


=  A  + 


ap 


Z(l-a)  * 


T ,13  +  a, 32  ° 


T ,14  4  a,42  0 


T , 22  ’  (T,2) 


1 , 23  =  ° 


T , 24  =  °  ’ 


4(l-a)  ’ 


-e2lV2x  “  2t  2.2  4  1T  2  ”  2^  4  2aT  22  4  2a  2T  2  =  A  4  \  P 


J 


where 


2  2 
v2  =  ^  +  32 


,  3,2  '  4,2 

(9x  )  (3x  ) 


At  this  point  it  seems  convenient  to  make  the  following 
substitutions 


-x 


> 


(6.16) 


(6.17) 


u  =  e 


<p  -  ua 


(6.18) 
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In  terms  of  the  new  variables  the  equation  becomes 


4(u-<{>) 

P  =  u  22 

u 


(6.19) 


Using  (6.19)  to  eliminate  p  from  the  ,  R  and  R0^  equations  we 

obtain  the  equivalent  system 


u,2*,2  "  u,lu,2  +  uu,12  ‘  uu,22 


l/r,2  .  2., 

-iV  x  +  u  A) 


(6.20) 


u  <J>  +  <j>  u  =  2uu  „  -  Au 

,22  >22  ,12 


(6.21) 


uu 


,11 


-  2  4)u 


,12 


u  22^u  “  2(^  +  \  v2“ 


+  u  -|  4*  9  ~  u  4)  =  0 


(6.22) 


By  use  of  compatibility  conditions  on  tj>  it  is  possible  to  obtain  a  third 
order  partial  differential  equation  involving  u.  However  it  is  so  unwieldy 
that  there  seems  to  be  little  hope  of  finding  a  general  solution.  Hence 
in  order  to  obtain  any  solutions  it  seems  reasonable  to  consider  ad  hoc 
simplifications.  Plausible  restrictions  might  be 


p  =  constant 


(6.23) 


p  =  constant  x(l-a) 


(6.24) 


u  -  u  =  0 


(6.25) 


Since  (6.25)  leads  to  the  greater  simplification  it  is  the  one  which  we 


shall  consider  here. 


i»i  <  tic  *  :  aud  - 
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From  equation  (6.25)  we  can  immediately  conclude  that 


t  1  ,  2  3  4, 

u  =  u(x  +  X  ,  X  ,  X  ) 


(6.26) 


However  we  also  know  that  to„=Tn/=0  which  implies  that 

,23  ,24 


,  3  4.  ,1,2. 

u  =  w(x  ,  X  )  v(x  +  X  ) 


(6.27) 


Returning  now  to  equation  (6.20)  we  can  write  $  in  terms  of  u  and  its 

2 

derivatives;  differentiating  with  respect  to  x  and  substituting  the 

resulting  expression  for  <J>  in  equation  (6.21)  we  obtain  as  an  explicit 

j  z  z 

formula  for  4>  in  terms  of  u, 


4> 


1  u  r„2  .  .  2n 

2  —2  [V  T  +  Au  1 

(u  ?) 

>  ^ 


(6.28) 


This  equation  gives  rise  to  a  new  expression  for  <f>  _  which  can  be  inserted 
in  equation  (6.20)  to  yield  the  differential  equation 


u  22^^t  +  Au2)  -  Au(u  ^  -  0  .  (6.29) 


Combining  (6.27)  and  (6.29)  we  get 


.  2  .  Av(v') 
Av  +  — 1 — L- 


v 


(6.30) 


The  left-hand-side  of  (6.30)  is  a  function  of  x  and  x  whereas  the 

1  2 

right-hand-side  is  a  function  of  x  +  x  .  Hence  it  follows  that  both 
sides  are  constant  and  we  get  the  two  equations 


V 


2 


T 


(6.31) 
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v" (Av2  +  k)  -  A v(v?)2  =  0  (6.32) 

where  k  is  a  constant.  To  solve  equation  (6.32)  let  v'  =  p,  v"  =  p  ^ 

Then  (6.32)  becomes  a  first  order  differential  equation  with  separable 
variables  and  its  general  solution  is 

v’  =  p  =  C(k  +  Av^)x/^  ,  (6.33) 

where  C  is  a  constant  of  integration.  Substituting  this  expression  for 
v'  in  equation  (6.32)  yields  the  equation 

v"  -  vC2A  =  0  (6.34) 


which  is  easily  integrated  to  give 

.  c/A  (x^  +  x2)  -C/K  (x^  +  x2)  (6.35) 

v  =  Ae  +  Be 

A  and  B  being  constants  of  integration.  Combining  equations  (6.33)  and 
(6.35)  we  obtain  the  following  expression  connecting  the  constants  k.  A, 

B  and  A  : 

k  =  -  4ABA  .  (6.36) 


Now  (6.27)  and  (6.28)  imply  that 

v  r n.2  ,  .  2  2n 

(j)  =  wv  -  - -  [V  t  +  Aw  v  ] 

2w(v'r 

Combining  this  result  with  (6.18),  (6.31)  and  (6.33)  we  get 


(6.37) 


(6.38) 
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In  order  to  solve  equation  (6.31)  note  that  u  = 

-T 

e  =  wv 

implies 


-  V2t  =  V2(log  w) 

(6.39) 

Introducing  the  substitution 


q 

w  =  e*  , 

(6.40) 

we  can  rewrite  (6.31)  in  the  form 


V2q  =  -  ke2q  . 

(6.41) 

Let  us  now  make  the  transformation 


3.4  -  3  .4 

z  =  X  +  IX  ,  z  =  x  -  IX 

(6.42) 

Substituting  in  (6.31)  we  obtain 


92q  k  2q 

9z9z  4  6 

(6.43) 

which  can  be  shown  [21]  to  have  the  general  solution 


2q  4  F' (z)  G' (z) 

k  [F(z)  +  G(z)J2 

(6.44) 

where  F  and  G  are  arbitrary  functions.  Combining  (6.40), 

(6.42)  and 

(6.44)  we  can  write 


2  4  F'(x3  +  ix4)  G'(x3  -  ix4) 

^  lr  L  Q  A9* 

K  [F(xJ  +  ix^)  +  G(XJ  -  ix  )] 

(6.45) 

Now  equations  (6.19)  and  (6.27)  imply  that 


■  .  ■ 


•  1  .  >  :  . 
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4 (1-a) v" 

P  =  - 

v 


(6.46) 


and  recalling  (6.34)  and  (6.38)  we  get 


p  =  2A 


(6.47) 


i.e.  in  the  universe  determined  by  this  solution  the  density  is  a  constant. 
We  have  thus  completely  specified  the  general  solution  of  equations  (6.16) 
when  they  are  subjected  to  (6.23).  For  completeness  we  express  the  results 
as  follows: 


-2x 


4  F'G' 
k(F+G) 


^c/ACx1  +  x2)  +  ^-C^x1  +  x2)1 


a  =  1 - -  ,  p  =  2A 

2C 


(6.48) 


J 


The  above  solution  is  valid  only  when  v?  and  V  t  are  not 
identically  zero.  If  we  now  assume  that  v'  is  zero  then  we  have  that 


3  ^ 

u  =  f  (x  ,  X  ) 


(6.49) 


Under  this  assumption  equations  (6.19)  -  (6.22)  become 


p  =  0 


„2  .  -2t 

V  t  =  -  Ae 


at/  3  4v 

4>  22  =  “  Af(x  »  x  ' 


(6.50) 

(6.51) 

(6.52) 


V  a  =  0 


(6.53) 


Integrating  (6.52)  we  get 


' 


■ 
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1  22  34  2-134  -1  34 

<f>  =  -  2  A(x  )  f  (x  ,  x  )  +  x  g (x  ,  x  ,  x  )  +  h(x  ,  x  ,  x  )  (6.54) 

which  implies  that 

1  2.2  ,  2  ,  1  3  4n  ,  ,  ,  1  3  4.  „  ,,N 

a  =  -  ^  A(x  )  +  x  g (x  ,  x  ,  x  )  +  h(x  ,  x  ,  x  )  .  (6.55) 

However,  a  32  =  a  42  =  ^  recluires  that  a  has  the  form 

a  ■  C(x'*',  x\  x^)  +  nCx1,  x2)  (6.56) 

and  comparing  (6.55)  and  (6.56)  we  conclude  that 

1  2.2  ,  2  ,  1*  ,  ,  ,  1  3  4N  ,, 

a  ■  -  ^  A(x  )  +  x  g(x  )  +  h(x  ,  x  ,  x  )  .  (6.57) 

Equation  (6.53)  implies 

V2h  =  0  .  (6.58) 


Equation  (6.51)  can  be  shown  to  have  the  general  solution 


s-2x  =  4  F1  (x3  +  ix4)  G*  (x3  -  ix4) 

5  A  3  4  3  4  2 

A[F(xJ  +  ix4)  +  G(xJ  -  ix4)] 


(6.59) 


F  and  G  being  arbitrary  functions  of  their  arguments.  Hence  the  complete 
solution  for  this  case  is: 


1  .  *  2 v  2  .  2,1.  ,  ,  ,  1  3  4. 

a  =  -  2  A(x  )  +  x  g (x  )  +  h(x  ,  x  ,  x  )  , 


A 


-2t  4  F'G' 


> 


(6.60) 


A  2  ’ 

(F+G)Z 


p  =  0 


J 


Let  us  now  assume  that  v'  is  not  zero  but  V  t  is  identically 
zero.  Equations  (6.31)  and  (6.32)  reduce  to: 


V2t  = 


0 


5 


(6.61) 
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w"  -  (v')2  =  0  ,  (6.62) 

where  it  has  been  assumed  that  A  ^  0.  From  equation  (6.61)  we  conclude 
that 

t  =  F(x2  +  ix^)  +  G(x2  -  ix^)  ,  (6.63) 

F  and  G  being  arbitrary  functions  of  their  arguments.  Integrating  (6.62) 
yields 

v*  =  av  ,  a  constant  ,  (6.64) 

which  in  turn  has  the  solution 

,  1  ,  2 

v  =  Aea  x  x  ,  A  constant.  (6.65) 

Returning  to  equation  (6.37)  we  can  show  that 

<p  =  (1 - -y)u  ,  (6.66) 

2a 


which  implies  that 


(6.67) 


Hence,  since  1  -  a  >  0,  we  must  have  that  A  >  0.  Defining  a  new  constant 
C  by 

a2  =  C2A  ,  (6.68) 


we  can  express  this  solution  in  the  form 

-2x  2C/A(x1  +  x2)  -  2F  -  2G 

e  =e 

a  =  1  -  — ,  p  =  2A  , 

2C 


_  (6.69) 

) 


. 


'  9 
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the  constant  A  having  been  absorbed  in  the  arbitrary  functions  F  and  G. 


Finally  if  both  V  x  and  v*  are  zero  then  the  equations  to  be 
satisfied  are: 


V  x  =  0  ,  p  =  0  ,  A  =  0  , 


\ 


a  n  =  a  „,  =  0  ,  a  no=  0  V  a  =  0 

,23  ,24  ,22  j 


(6.70) 


The  general  solution  to  this  set  is 


-2x  -2(F+G)  n 

e  =  e  ,  p  =  0  , 


2  (  1.  ,  ,  ,  1  3  4. 

a  =  x  g (x  )  +  h(x  ,  x  ,  x  )  , 


(6.71) 


where  g  is  an  arbitrary  function  of  x  and  h  is  an  arbitrary 

3  4 

harmonic  with  respect  to  the  variables  x  and  x  . 


function 
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CHAPTER  VII 

Examples  and  Discussion  of  Results. 

Let  us  now  investigate  some  particular  solutions.  Turning  to 
(6.48)  we  observe  that  the  two  arbitrary  functions  F  and  G  may  be  chosen 


to  simplify 

the  form  of  the  solution.  One  possibility  is  to  put 

3  4  3  4 

F  =  x  +  ix  -  R  ,  G  =  x  -  ix  -  R  ,  (7.1) 

where  R  is  a  constant.  The  arbitrary  constants  A  and  B  may  be  chosen 


such  that 

A  =  B  =  |  ,  (7.2) 

and  substituting  in  (6.36)  we  get 

k  =  -  A  .  (7.3) 


The  general 

solution  (6.48)  thus  reduces  to 

-2t  1  ,  2„  rr,  1  .  2X 

e  -  „  „  cosh  C/A(x  +  x  ) 

A(xJ  -  Rj 

(7.4) 

a  =  1  -  — ,  p  =  2A  . 

2C 

Applying  the  coordinate  transformation  (6.3)  and  substituting  in  (6.4) 


we  obtain  the  metric  tensor 


' 


r-n  .^*1  't  .1  c<  •  n  It  is«; 
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gij  = 


r 


2C‘ 


-  1  + 


2C 


0 


0 


V 


-  i  + 


2  - 


2C‘ 


2C 


0 


0 


0 


0 


Cosh^c/2/T  x2 

-3  2 

A(xJ  -  R) 


0 


0 


A 


0 


0 


Cosh2Cv/2A  x2 

A(x3  -  R)2 


y 


(7.5) 


At  the  origin  of  the  coordinate  system  it  would  be  desirable  if  (7.5)  at 
least  approximates  the  Minkowski  metric.  This  would  be  true  if  we  put 


1 

2  ’ 


"  “  R2 


(7.6) 


However  in  discussing  the  Einstein  static  universe  we  observed  that  A  had 
the  order  of  magnitude  of  the  inverse  square  of  the  radius  of  the  universe 
and  we  here  obtain  the  same  result  if  R  is  interpreted  as  the  radius.  If 
the  relations  (7.6)  are  assumed  to  be  exact  the  metric  tensor  becomes 


gij  = 


/ 


-  1 


0 


0 


0 


0  0 


0 


0 


\ 


Cosh2 (x2/R) 

4  - 


0  0 


0 


0 


0 


Cosh2(x2/R) 

-3  0 

<T-X)  J 


(7.7) 
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Identifying  x  with  a  radial  coordinate  it  seems  reasonable  to  regard 
(7.7)  as  characterizing  a  static  universe  with  a  constant  mass  density 
and  with  a  spatially  closed  geometrical  structure.  Note  that  for  this 
model  A  must  be  non-negative. 

The  solution  (6.69)  yields  a  particularly  simple  example  when 

we  choose 


F  =  G  =  0 


(7.8) 


We  then  obtain 


a 


-2x  2C/A(x^  +  x2) 

e  =  e 

±  1 
2C2 


2A  . 


As  in  the  previous  case  it  is  desirable  to  choose  C 
may  be  written 


>  (7.9) 

y 

The  metric 


giJJ 


-  1 

0 

0 

0 

V 


0  0 


0 


\ 


10  0 
0  e2^  X  0 


0 


0  e 


2/A  x2 


/ 


(7.10) 


We  observe  that  in  this  universe  the  only  singularity  is  at  infinity  so 
that  it  would  seem  to  correspond  to  a  spatially  infinite  universe  with 
constant  density. 


g(x1) 


=  0 


Finally  let  us  consider  (6.71)  where  we  put 

3 

,  h  =  h(x  )  so  that  the  solution  becomes 


F 


G  =  0, 


' 
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-2t 

e 


1 


a  =  h(x^) 


p  =  0 


(7.11) 


and  the  line  element  may  be  written  in  the  form 

ds^  =  [-  1  +  h(x^)](dx^)2  -  2h(x^)dx^dx^  +  [1  +  h(x^)](dx^)^ 

-32  -42 

+  (dxJ)  +  (dx4) 


(7.12) 


If  we  now  identify  x  with  <j>  ,  x  with  r  and  x  with  z,  where  r,  <j) 
and  z  are  cylindrical  polar  coordinates,  then  there  is  a  similarity 
between  the  line  elements  (7.12)  and  (5.14).  Recall  that  the  latter 
describes  a  flat  space  rotating  about  the  axis  of  z  with  constant  angular 
velocity  oo  .  From  (6.6)  it  is  easy  to  obtain 


u.  =  -  (1  -  h  ,  h  ,  0  ,  0) 

1 


(7.13) 


where  u^  is  the  covariant  velocity  4-vector.  Straightforward  computation 
yields 


( 

k  3 

- ?  '/9'  for  i  =  1  ,  j  =  3  , 

2(1  -  h)X/ 

h  A2  -  h) 

-  — * - —rr  for  i  =  2  ,  j  =  3  ,  (7.14) 

2(1  -  h)  ' 

0  otherwise, 

V 


and  it  can  further  be  shown  that 
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/ 


u 


j  ,k] 


0 


even  permutation  of  1,  2,  3  , 


odd  permutation  of  1,  2,  3  ,  (7.15) 


otherwise 


V 

For  the  metric  (7.12)  the  vector  ft1  defined  by  (5.15)  can  be  readily 
computed  from  (7.15)  to  give 


(0,  0,  0,  - 


(7.16) 


By  analogy  with  (5.18)  the  angular  velocity  of  the  bulk  matter  in  this 
universe  is 


h 


0)  =  - 


2(1  -  h)  * 


(7.17) 


-3  2 

However  combining  h  =  h(x  )  and  V  h  =  0  we  find  h  must  have  the  form 


h  =  ax  +  b 


(7.18) 


where  a  and  b  are  constants.  If  we  require  h  =  0  at  the  origin  we  get 
b  =  0  and  we  may  write 


-3 

,  x 

h  "  R 


(7.19) 


where  R  =  -  .  Substituting  (7.19)  in  (7.17)  we  obtain 

cl 


CO  = 


2(x3  -  R) 


(7.20) 


Hence  for  this  model  the  angular  velocity  depends  on  radial  distance  from 
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the  origin  though  it  is  independent  of  the  timelike  coordinate  x 
From  (7.20)  it  is  evident  that  there  is  a  singularity  at  r  =  R. 


Combining  (7.12)  and  (7.19)  we  may  write  the  line  element  in 

the  form 


ds2  =  [-  1  +  —]  (dx1)2  -  2  “  dx1dx2  +  [1  +  ~](dx2)2 

-32  -42 

+  (dx"5)  +  (dx  ) 


(7.21) 


Let  us  now  investigate  the  behaviour  of  light  in  the  universe 
described  by  (7.21).  We  recall  that  light  rays  travel  along  null  geodesics 
for  which  ds  =  0  and  we  may  thus  write,  for  the  velocity  of  light  travelling 
radially  towards  a  distant  nebula. 


dr  f  r.  1/2 

dt  U  R} 


(7.22) 


-1  -3 

where  for  convenience  we  write  t  for  x  and  r  for  x  .  From  (7.22)  we  get 


e  r 


t  = 


dr 


J0  (1  -  H)1/2 


(7.23) 


Hence  if  t  is  the  time  of  emission  of  light  from  a  nebula  located  at  the 
radial  distance  r  and  t^  is  the  time  of  reception  by  an  observer  at  the 
origin  we  obtain 


t2  =  tl  + 


dr 


0  (1  -  |)1/2 


(7.24) 


from  which  we  deduce  that 


At2  =  [1  + 


(i  -  g) 


1  dri  A  M 

r.1/2  dt  Atl 


9 


(7.25) 


■  L  1  ,-jM 
. 
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where  At^  is  the  period  of  the  emitted  light  wave,  At^  is  the  period  of 

dr 

the  received  wave  and  —  is  the  radial  velocity  of  the  nebula  at  the  time 
of  emission.  However  the  proper  time  intervals  (At-)  and  (At,,)  are 

1  p  2  p 

given  by 


(Atpp  -  (1  -  |)1/2  Atx 


(7.26) 


(At2)p  =  At2 


(7.27) 


Hence 


A  +  AA 


<Atl>p 


(1  - 


U-|) 


dr 

dt 


(7.28) 


and  the  expression  for  wavelength  shift  is 


z  = 


A A  _  1 _  1  dr 

x  ”u  - 1)1/2  u-|)  dt 


-  1 


(7.29) 


Suppose  now  that  the  light  source  has  zero  initial  radial  velocity  with 
respect  to  the  observer  at  the  origin.  The  expression  (7.29)  becomes 

z  =  (1  -  '  1/2  -  1  .  (7.30) 

If  we  assume  that  the  magnitude  of  R  is  very  large  then  a  Maclaurin's 

r  —  ]_/  2 

series  expansion  of  (1  -  -)  ,  neglecting  second  order  terms, 

R 

yields 

z  =  a‘r  ■  (7-31) 

so  that  in  first  approximation  the  velocity  distance  relation  is  satisfied. 


as  long  as  r  <<  R  . 


. 
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The  line  element  (7.21)  thus  appears  to  describe  a  universe  which 

-4 

rotates  about  the  x  -axis  and  which  also  gives  rise  to  a  wavelength  shift 
in  light  rays  which  travel  radially  from  a  distant  nebula  to  the  axis  of 
rotation.  We  note  that,  like  the  de  Sitter  model,  this  universe  is  empty 
and  the  nebulae  must  be  regarded  as  test  particles  of  negligible  mass.  It 
is  to  be  hoped  that  more  sophisticated  solutions  may  be  found  which  will 
not  impose  such  severe  restrictions  on  x  and  which  will  retain  the  expanding- 
rotating  character  of  (7.21)  while  admitting  a  positive  mass  density. 
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